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SUMflARY 


An algorithm is described for selecting a grid subset for calcu- 
lating radiative transport. The subset spacing is determined by using 
the variation in aerothermai properties across the full grid of the 
shock layer. Results show that a radiation grid subset of 15 to 20 
points can be used for viscous-shock-layer calculations where approx- 
imately 50 grid points are required to define the aerothermal profiles. 
Results are presented for various planetary entry conditions with and 
without mass injection to lemonstrate both the validity and utility 
of the algorithm. 




' - . ■■— - ^ -I 

/j- 

i 

. • ! 
<L 

r 


SYMBOLS 
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p. _ 1 
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Cj 

mass fraction of species j 
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h t 

*“ 1 
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C' 

sum of the absolute values of the species concentration mass 
fraction normal derivatives 
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e' 

radiative flux divergence (normal derivative) 


m 

number of radiation grid subset points 
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nondimensional mass injection rate, (pV)^(pV)^ 
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Np.e 

Reynolds number * 
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number of chemical species 
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n 

distance from wall or number of shock-layer grid points 
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shock standoff distance 

temperature »- 
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temperature normal derivative 
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weighting factor 
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quantity defined by equation 1 
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species concentration 
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INTRODUCTION 


Prediction of wall heating and injection rates associated with a hypersonic 
planetary entry vehicle requires a detailed analysis of thermodynamic and trans- 
port processes in the high energy flow field surrounding the vehicle. The 
most costly part of the analysis is typically the calculation of radiative 
transport, which may consume 80 to 90 percent of the total computer solution 
time. Since viscous-shock-layer solution techniques, such as those of Moss 
(ref. 1), and Anderson and Moss (ref. 2), require 50 to 100 grid points across 
the shock layer to adequately resolve all significant transport events, com- 
putation of radiative transport at all grid points would make these techniques 
prohibitively expensive for routine usage. 

The present report describes an algorithm which uses the variation in 
aerothermal properties across the shock layer to select a grid subset for 
computing the radiation transport. Test cases run to validate the algorithm 
are discussed. Results are presented for Earth and Jupiter entry conditions 
with moderate to massive injection. To indicate the stability of the algorithm, 
results are presented for cases in which the full grid does not adequately 
model important transport events. 
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PROBLEM 


Accurate computation of radiative transport within a high temperature 
shock layer, such as that experienced by a planetary probe during entry, requires 
detailed consideration of both continuum and line transitions for ionic, atomic, 
and molecular species within the layer. Because the strengths of the various 
transitions are functions of temperature and species concentrations, strong 
nonlinear couplings can develop between the radiation transport processes, the 
other shock -layer transport processes, and the mass injection at the wall. 

Figures 1 and_2 show the flux divergence and temperature profiles for 
typical entry cases and illustrate some of the difficulties encountered in 

4 - 

selecting a radiation transport grid for the shock layer. The shock layer- 
consists of a weakly viscous outer region of atmospheric species that have 
been dissociated and partly ionized by the bow shock to form a strongly 
emitting plasma, a relatively cool injection region near the wall, and a strong- 
ly absorbing intermediate viscous mixing layer where the injection. species are 
dissociated by the shock energy. A successful radiation subgrid selection 
algorithm must identify and properly weight complex events occuring over very 
small segments of the shock layer. 
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METHOD 


Since radiation transport within the shock layer is dependent primarily 
upon the temperature and distribution of the constituent species, intuition 
suggests that the rates of change qf these quantities, normal to the wall, 
should provide sufficient information to define a suitable radiation grid. 
Numerical experiments indicate that the absolute values of the temperature 
derivatives and the sum of tbe absolute values of the species concentration 
derivatives can be used to identify the spatial regions where the major radia- 
tion transport events occur. This is illustrated in figures 1 and 2. For 
these cases, dissociated injection species form a strongly absorbing region 
near the wall, which coincides with a peak in the concentration derivative. 
Large changes in the divergence of the radiation flux within the mixing region 
are also marked by the concentration derivative. The temperature derivative 
provides additional information about the mixing layer and its decay marks 
the transition to the emitting outer layer. 

The radiation subgrid is therefore defined by weighting each point of 
the full shock-layer grid according to the derivatives of temperature, concen- 
tration, and, if applicable, the flux divergence computed at the previous 
iteration. (Including the divergence damps out possible oscillations in the 
subgrid point selection procedure during the first few iterations of the 
flow-field solution.) The grid subset is chosen In equal increments of Z, as 
follows. If 
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I = L 2 ^ m) (2) 
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The method is illustrated in figure 3. To avoid interpolation when 
Z|^<Zji< Zk+i the subgrid point is chosen to coincide with the point 
k of the full grid. 

The spacing of the radiation subgrid points is controlled by the mag- 
nitudes of the weighting factors W. If the weighting factors are set to 
zero, then Z = n, and every (g) point will be selected. (For ej<ample, 
if n = 50, m = 10, every fifth point will be selected.) A high value of 
Wj. will tend to concentrate grid points in spatial regions where the 
species concentrations are changing rapidly. A high value of will con- 
centrate grid points in the mixing region. A high value of We will cluster 
grid points wherever the flux divergence is rapidly changing. 
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RESULTS AND DISCUSSION 


The radiation grid subset selection procedure defined by equations 
1 and 2 and the effects of the weighting factors of equation 1 were 
evaluated using the Earth, Venus, and Jupiter entry environments described 
in Table I. Cases run included no injection, atmospheric element injection, 
coupled ablation* and prescribed ablation injection. For each condition, 
a benchmark -case stagnation profile was obtained for fifty points across 
the shock layer. The hypersonic viscous-shock-layer program (HYVIS), 
which is based on the analysis of references 1 and 2, was used for all 
computations. From the resulting chemistry and temperature profiles, the 
radiative flux divergence was recomputed for grid subsets of 20, 15 and 
10 points and compared with the benchmark results. For each injection 
condition, at least one case was run in which the complete stagnation 
profile was^ computed using a radiation grid subset. 

Two radiation models were used. Initial testing focused upon Earth 
entry using the Engel radiation code (ref. 3) which uses the coupled line- 
continuum model of Wilson (ref. 4) to compute the flux divergence directly. 
Using a 20 point radiation grid and weighting factors of unity, stagnation 
profiles were computed for cast: of no injection, air injection, and 
coupled ablation product injection. A local quadratic interpolation 
method was used to estimate the flux divergence over the full shock- layer 
grid. For ablation cases, all state and transport properties agreed with 
benchmark solutions to within 2.5 percent. Cases with no injection and 
air injection agreed to within 1.0 percent. 
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Further testing was conducted with the Aerotherm -69 radiation code 
of Nicolet (ref. 5) which computes the bidirectional radiative flux at 
each grid point. The HYVIS program computes the divergence by differen- 
tiating the net flux profile using a three-point differencing technique. 
Initial results with the Aerotherm code showed unstable behavior of the 
flux impinging upon the wall and the flux divergence at both the wall and 
the shock as the number of points in the radiation grid was reduced. To 
alleviate this instability, spline interpolation and differentiation for 
the radiative flux and divergence were tried, and the weighting factors 
were varied, all without success. The instability was removed by 
forcing the radiation subgrid to include the two shock layer grid points 
nearest the wall and two at the shock to “tie down" the ends of the flux 
divergence curve. 

Table II shows results obtained by recomputing the flux divergence 
for the case of figure 1 with radiation subgrids of 20, 15, and 10 points. 
Four methods are compared. The basic method used equations 1 and 2 to 
pick the grid subset and quadratic interpolation and three point dif- 
ferencing to compute the flux divergence. In the second method, spline 
interpolation and differencing were used. The third method tied down the 
four end points but used interior subgrid points hand selected by an 
analyst to represent overall property variation within the layer. (This 
is equivalent to setting the weighting factors to zero.) For the fourth 
method the interior subgrid was randomly selected. 

For this moderate blowing case, the first three procedures yielded 
satisfactory results for subgrids down to 15 points. Using 10 points, 
none of the methods was completely satisfactory. The spline technique 



perfomed marginally ''•-'■;ter in the interior of the layer but did poorly 
at the end points. This and other cases indicated that use of the local 
quadratic interpolation scheme with three -point differencing is less 
likely to cause gross errors if the selection of the sub grid does not 
adequately model the radiation transport events. 

The behavior of the basic method when varying the weighting factors was 
also investigated. Weighting factors less than 0.5 shifted the subgrid too 
much towards the shock, thus losing detail in the mixing region. Weighting 
factors greater than 1.5 concentrated poi*’ts too heavily around extremum 
points of the flux divergence profile. Weighting factors of unity were 
found to be optimum for the cases considered. 

Table III shows stagnation profile data for Jovian entry with pre- 
scribed massive injection. Cases were run for grids which were both 
course (an = 0.02ns) and fine (a n = 0.005ns) in the region adjacent to 
the wall. For each grid a twenty-point radiation subgrid was used to 
obtain an estimate of the stagnation profile. Then this estimate was 
used to start the fifty point solution. 

The wall radiative heating was adequately predicted for all the cases, 
the total spread being less than 5 percent. Considerable discrepancies 
appear, however, in the flux divergence and the nonradiative (conduction 
plus convection plus diffusion) heating results. Examination of the 
temperature and species concentration profiles indicated that the coarse 
wall grid spacing did not adequately model the absorption layer even 
with a full radiation grid. With the more appropriate fine wall grid 
spacing, a twenty point radiation subgrid adequately modeled the wall 
radiative heating and flux divergence. The wall nonradiative heating 
also showed excellent agreement. 
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Cases were also run for Venus entry, with no injection and moderate 
injection^ using radiation subgrids of 20, 15, and 10 points. The trends 
were very similar to those for earth entry (Table II), so the results are 
not presented. 
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CONCLUSIONS 


An algorithm is presented for automatically selecting a radiation 
grid subset based on the variation of aerothermal properties across a shock 
layer. The method has been tested with two radiation models, three 
planetary entry conditions, and a wide range of blowing rates. Results 
presented show that the method yields accurate predictions of radiative 
transport over a wide range of injection conditions with 15 to 20 
subgrid radiation points provided the full grid (50 points), adequately 
models the major transport processes. Since the cost of oerforming the 
radiative calculation is almost directly proportional to the number of 
grid points, the present algorithm provides a means for substantially 
reducing the computational cost for calculating radiative flow fields. 
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Altitude, km I 

Velocity, km/sec 
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Mach number 
Temperature, k 
Density, g/cm^ 

Nose radius, m 

Atmospheric Composition, 

Mass Fractions 

Injectant Composition, 

Mass Fractions 
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ENTRY CONDITTONS FOR 


lATION SUBGRID TEST CASES 


Earth 
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TABLE I I. -COMPARISON OF RADIATION SUBGRID MODELS 


Flux at 
Wall 
MW/m2 


Divergence 
at Wall, 
GW/m3 


Divergence 

Minimum, 

GW/m3 


Di vergence 
at Shock, 
GW/m3 


Baseline 
(50 points) 

8.43 

-1.259 

i 

! 

-1.707 

23.47 

Basic Method: 





20 points 

8.26 

-1.239 

-1.679 

23.54 

15 points 

8.32 

-1.254 

-1.668 

23.37 

10 points 

8.25 

-1.336 

-1.534 

22.72 

Splines: 





50 points 

8.43 

-1.223 

-1.720 

21.40 

20 points 

8.23 

-1.215 

-1.693 

20.30 

15 points 

8.11 

-1.227 

-1.710 

19.83 

10 points 

7.91 

-1.210 

-1.656 

19.98 

Hand Selected 





20 points 

8.25 

-1.254 

-1.698 

23.49 

15 points 

8.30 

-1.264 

-1.686 

23.37 

10 points 

8.32 

-1.327 

-1.482 

11.75 

Random Grid 





20 points 

8.24 

-1.223 

-1.576 

23.76 

15 points 1 

8.13 

-1.191 

-1.524 

23.79 

10 points 

7.75 

-1.037 

-0.991 

24.33 
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TABLE III. -JOVIAN ENTRY WITH PRESCRIBED INJECTION 
(Ni>e = 1.34x10®, p = 9.86 atm, m = 0.40) 


CW6-50 


Grid Points 
Wall Mesh 
Shock Standoff, m 


Wall radiative 

heating, Mw/m2 

Continuum 

Line 

Total 


Flux Divergences, 

TW/m3 

Wall 

0.05*ns 

0.88*ns 

1.00*ns 


Wall nonradiative 
heating, MW/m 


0.01896 


191.4 

58.8 

250.1 


20 

Coarse 

0;01719 


188.8 

61.0 

249.8 


-2.320 


-4.073 


FW6-50 

FWG-20 

50 

20 

fine 

fine 

0.01877 

0.01390 

194.2 

192.3 

64.7 

69.9 

258.9 

262.2 

-65 

I -62 

-290 

-291 

116 

113 

503 

500 

0.908 . 

0.849 
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Figure 1. - Radiative flux divergence and other selected 
shock-layer properties (Earth entry typical) . 
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* To avoid interpolation, move 
S. to coincide with k. 



Figure 3. -< Radiation subgrid selection procedure 
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